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A deep X-ray observation of M82 with XMM- Newton 
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ABSTRACT 

We report on the analysis of a deep (100 ks) observation of the starburst galaxy M82 
with the EPIC and RGS instruments on board the X-ray telescope XMM-Newton. 
The broad-band (0.5-10 keV) emission is due to at least three spectral components: i) 
continuum emission from point sources; ii) thermal plasma emission from hot gas; iii) 
charge exchange emission from neutral metals (Mg and Si). The plasma emission has a 
double-peaked differential emission measure, with the peaks at ~ 0.5 keV and ~ 7 keV. 
Spatially resolved spectroscopy has shown that the chemical absolute abundances are 
not uniformly distributed in the outflow, but are larger in the outskirts and smaller 
close to the galaxy centre. The abundance ratios also show spatial variations. The 
X-ray derived Oxygen abundance is lower than that measured in the atmospheres of 
red supergiant stars, leading to the hypothesis that a significant fraction of Oxygen 
ions have already cooled off and no longer emit at energies > 0.5 keV. 

Key words: galaxies: individual: M82 - galaxies: abundances - X-rays: ISM - ISM: 
jets and outflows - plasmas - atomic processes. 



1 INTRODUCTION 

M82 is a nearby galaxy 

type starbu rst galaxy (iRieke et alj Il980t ). Its proximity 
(3.63 Mpc, iFreedman et alj Il994h and luminosity (10 44 

erg s" 1 in t he far infrare d, 10 40 erg s _1 in the X- 

ray domain; iRanalli et alj I2003T ). along with the pres- 
ence of an ongoing powerful (~ 3 Mq yr _1 ) starburst 
has made it the subject of countless studies at all wave- 
lengths. The starburst, located in the galaxy central re- 
gions, is driving a several kp c large outflow perp endicu- 
lar to the plane of the galaxy l|Heckman et al.lll990l ) which 
has been observed at multiple wavelengt hs, such as the ra- 
dio dSeaquist fc qdegardlll99ll ). infrared (| Alton et all 1 19991: 



Engelbracht et al.1 I2006D , Ha and X-rays ( Lehnert et al.l 
19991 ). Recent observations in the X-ray domain with the 



Chandra, XM M-Newton a nd Su z aku observatories are re - 
ported in iGriffiths etaH d2000h. iRead fc Stevens! (|2002h . 



lOriglia et ahl (|2004h , IStrickland fc Heckmanl (|2007l ) and 
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ITsuru et al.l l|2007l ). The small inclination of the galaxy disc 
makes M82 observed almost exactly edge-on, thus allowing a 
good perspective on the supernova-driven outflow and con- 
tributing to its success as an astrophysical laboratory. 

M82 was thus chosen as a test case where the metal 
enrichment due to a burst of star formation can be di- 
rectly measured. Determining the abundance of key ele- 
ments released in the interstellar medium (ISM) by super- 
novae (SN) with different mass progenitors and hence on 
different time scales, will have a strong astrophysical im- 
pact in dr awing a global p i cture of galaxy format ion and 
evolution (|McWiniamlll997l ; iMaeder fc Contilll994r ). While 
the metals locked into the stellar atmospheres give a pic- 
ture of the abundances at the beginning of the last burst of 
star formation (SF), the hot gas heated by core-collapse su- 
pernova explosions and emitting in the X-rays should trace 
the enrichment by the SN which were formed in the new 
generation of stars. 

In our previous paper l|Origlia et alj|2004h we presented 
high resolution infrared (J and H band) spectra that we ac- 
quired with the 3.6 m Italian Telescopio Nazionale Galileo 
(TNG), together with a rather shallow (~ 15 ks) X-ray ob- 
servation obtained from the XMM-Newton archive. Our re- 
sults hinted for a confirmation of the expected scenario in 
which the gaseous component has a higher content of a- 
elements than the stellar one, and a similar content of Fe. 
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However, some new issues were posed, since we found a very 
low abundance of O and Ne with respect to other a-elements 
(e.g., O/Mg ~ 0.2, Ne/Mg ~ 0.3) in the hot gas present in 
the central (< 1 kpc) regions of M82, which could not be 
satisfactorily explained. For these reasons we were granted 
the deeper (~ 100 ks) observation of M82 which we present 
in this article. Two re ports on an early a nalys is appeared 
in IRanalli et all d2005l) and IRanalli et all (|2006l ). Both this 
work and lOriglia et al.l (|2004T ) are part of an ongoing larger 
effort to measure metallicity enhancements in a sample of 
galaxies. 

The X-ray spectrum of M82 is rich and complex. 
For a review about the physical and observational pa- 
rameters which may be recovered through X-ray spec- 
troscopy (shape of the spectra; slopes of power-laws; tem- 
peratures and abu ndances of thermal pla smas; etc.), we 
refer the reader to iPaerels fc Kahnl |2003), while for the 
physical processes of thermal and charge- e xchan ge emis- 
sion we refer to th e book by Ostcrbrock (|l989l) and to 
the recent review bv lKallman fc Palmeril ( 20071 ); for plasma 
emission we mention both histor ical references such as 
lLandini fc Monsignori Fossil (|l970l ) and the current ones 
whi ch illustrate the st ate of the art fe.g.. iMewe et af1ll995l 
and lSmith et al.ll200ll ). 

This paper is structured as follows: Sect. [2] presents the 
XMM- Newton observation and data reduction of M82. In 
Sect.[3]the EPIC spectra of the central region are analysed, 
and results are shown about the temperature distribution 
of the plasma, charge exchange emission, and a comparison 
of different plasma codes. In Sect. [3] the RGS spectra are 
analysed, improving and extending the results obtained from 
EPIC data. In Sect.[5]a spatially-resolved spectral analysis of 
the EPIC data of the M82 outflow is presented; a dependence 
of the chemical abundances on the height from the galaxy 
plane is found. The implications of the main findings of this 
article are discussed in Sect.JS] Finally, in Sect.[7J the results 
are summarised. 

Throughout this paper, the abundances a re linearly 
scaled according to the lGrevesse fc Sauva 3 (|l998l ) solar com- 
position (i.e., we always mean X/Xq for each chemical ele- 
ment X, following the common usage of X-ray based litera- 
ture). We assume 3.63 Mpc as the distance to M82 , based 
on Cepheids observations by Frccdma n et al.l l| 19941 ). 



2 OBSERVATIONS 

M82 was observed with XMM-Newton on April 21 st 2004 for 
about 100 ks (Obs-id: 0206080101). After screening for back- 
ground flares, about 73 ks of data were accepted for analy- 
sis. A true-colour image of M82 is showed in Fig. [T] with the 
red, green and blue colour channels representing the 0.4- 
1.0 keV, 1.0-2.0 keV and 2.0-8.0 keV bands, respectively. 
The images in the three bands have been smoothed with 
the asmooth SAS task , and scaled for the best visual pre- 
sentation according to iLupton et al.l (|2004l ). Superimposed 
on the image, we also show the regions from which the EPIC 
spectra were extracted. The RGS spectrum covers an area 
centred on the galaxy nucleus and extends approximately 
over the region shown in Fig. [1] as the blue dashed parallel 
lines. 



SPECTROSCOPY OF THE CENTRAL 
REGIONS OF M82: EPIC DATA 



Th e EPIC data for bot h the MOS (|Turner et al.ll200lh and 
pn l|Striider et al.ll200ll ) cameras were analysed with the SAS 
(version 6.5.0) and XSPEC (version 11.3.2z) software. We 
began the spectral analysis from the central region of M82, 
by extracting spectra from a circle with diameter 1' and cen- 
tred on the coordinates 09:55:51 and +69:40:39. In this area 
the largest emission from hot plasma is pr esent; however, 
many point sources exist in the same region ([Griffiths et al.l 
2000). Thus, the spectral model should account for both 
plasma emission and point sources, along with considerable 
absorption. 



3.1 Spectrum of underlying point sources 

To define an optimal model for the point sources, we anal- 
ysed the spectrum in the 3-8 keV band, and found that it can 
be described by a power law with photon index T = 1-601q'o3 
(all errors on spectral parameters are at 90%), and X 2 =682 
(with 571 degrees of freedom, hence %rcd ~ 1-19). This emis- 
sio n may be ascribed, for all practical purposes, to M82 X- 
1 (jMatsumoto et all 1200 ll ) because the lu minosity of this 
sourc e is much larger than the other ones l|Griffiths et al.l 
2000). This best-fitting value for V was not kept frozen in 
the spectral fits discussed below, but rather used as a con- 
sistency check. The point sources, although accounted for 
in the spectral fits, will not be discussed further. Results 
ab out M82 X-l from this observation have been presented 
in iMucciarelli et all (120061 ). M82 X-l is also clearly visible 
in Fig. \T\ as the blue spot near the galaxy centre. 



3.2 Background spectra for extended sources 

The background spectra were derived from blank-sky obser- 
vations. Because of the source extent, and the spatial de- 
pendence of the background spectrum on the XMM-TVewfon 
CCDs, it is in fact impossible to obtain background fields 
from the M82 observation. While background spectra could 
be extracted near the edges of the field of view where the 
M82 emission is negligible, these spectra would be different 
from what is expected in the detector region covered by M82. 
Thus we used the XMM - Newton blank-sky files provided by 
the Birmingham group (Re ad fc Ponmanll2003l ). 

Blank-sky observations have the advantage that a back- 
ground spectrum can be extracted in the same detector re- 
gion covered by the source, but they need greater care in 
checking that the background level be as close as possible to 
the source's one. In order to match the background bright- 
ness, we selected two large circular regions (radius ~ 2.35') 
in a position several arcmin off-axis where no emission from 
M82 is present, and extracted spectra from both the M82 
data and blank-sky files. In this way, spectra extracted from 
the same detector region are compared in order to scale 
the brightness of the latter to the former's. We found only 
small differences: the blank-sky /source flux ratios in these 
two regions are 0.97 ± 0.03 (MOS1), 0.91±g g| (MOS2) and 
0.86^Iq'q2 (P n )i an d are largely independent from the details 
of model fitting. Thus, in the following analysis we rescaled 
every background spectrum by means of the BACKSCAL 
keyword in the FITS files. 
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Please insert fig01.jpg here. 

Figure 1. Spectral regions superimposed on XMM- Newton image. The cyan boxes correspond to the outflow regions analysed in Sec. [5] 
while the small yellow circles show the areas around the point sources which have been excluded from the extraction of the spectra. The 
large white circle shows the central region whose EPIC spectrum is discussed in Sec. [3] because of the large number of point sources 
found in this region (not evident in the image, and not indicated by circles to avoid cluttering the picture) the area around them has 
not been excluded, but rather the point sources have been accounted for in the spectral model. The blue dashed lines sketch the RGS 
spectral extraction region. They are parallel to the dispersion direction, and are spaced by 2' which corresponds roughly to the 90% 
of the PSF width, the value used in Sect. [4] for the spectral extraction. Along the dispersion direction, the PSF is only limited by the 
instrument vignetting. 



3.3 Temperature distribution of the plasma 

It is customary to model the gaseous emission in star form- 
ing galaxies with one or more single-temperature plasmas. 
A ' warm' component around 0.7-1 keV is often reported 
(e.g. iPtak et all 1 19991 ). A 'hot', high en ergy component is 
also sometimes present ijdella Ceca et al1ll997l : ICappi et all 
1 19991 ) in studies based on ASCA or BeppoSAX data, de- 
scribed either as a power-law, or as a thermal spectrum 
with kT ~ 4-10 keV. Because of the poor spatial resolu- 
tion of those satellites it is not clear if this component is re- 
ally diffuse or rather due to point sources. Thus, some care 
should be taken in comparing literature results to this paper, 
since in the following we will account for both point sources 
(via the power law discussed in the previous paragraph) 
and for hot gas. A 'cold' componen t, around 0.1-0.3 keV , 
has been also occa sionally observed (|della Ceca et al"lll999l ; 
iTsuru et all l2007f ). Moreover, a po tential pitfall has bee n 
discovered (the so-called 'Fe-bias'. iBuote fc Fabianl Il998h . 
which occurs when fitting single-temperature models to data 
which are intrinsically multi-temperature, and which leads 
to artificially low abundances. 

The use of multi-temperature models that allow a fine 
control of the Differential Emission Measure (DEM) and are 
also insensitive to the Fe-bias might be preferable. However, 
this method has been seldom used because of its demanding 
requirements in terms of data quality and computing time. 
The quality of our M82 observation is sufficiently high to 
allow for the first time the use of multi-temperature models 
even in spatially-resolved spectroscopy; 

Tw o spectral codes, M EKAL (|Mewe et all Il995l ) and 
APEC (|Smith et all l200ll ). are currently used in the lit- 
erature. However, while several versions of the common 
MEKAL code are present with multi- temperature flavours 
jLemen et all Il989l ; ISingh et all Il996l ) in the XSPEC dis- 
tribution, no multi-temperature APEC model is present. 
Since we intended to test which code offers the best per- 
formances, we built multi-temperature APEC models for 
XSPEC, which behave like the MEKAL ones, with the only 
difference that the APEC routines are called. 

We used initially the XSPEC c6pvmkl model, which 
parametrizes the plasma DEM (i.e. the amount of plasma 



DEM sampling was already sufficiently good: as it may be 
seen, the high-temperature component is distributed like a 
rather narrow bell curve, so little difference is found between 
a single temperature element and the sum of a few elements 
with very similar temperatures. 



3.4 Results and comparison of plasma codes 

The abundances for O, Ne, Mg, Si, S, Fe were left free to 
vary, since these elements have strong lines in the considered 
wavelength range. Other a-elements with weaker lines (Na, 
Al, Ar, Ca), which would not be sufficiently constrained on 
their own, were initially left at the Solar value. However, 
after noticing that the metals close in atomic number have 
similar abundances relative to the solar values (Table |TJ), we 
hooked Na and Al to Mg, and coupled Ar and Ca together 
in order to achieve a better x 2 - An F-test showed that this 
approach, as opposed to leaving Na, Al, Ar and Ca to solar 
values, leads to an improvement of the x 2 with significance 
> 99.5%. While this can be regarded as low significance, we 
still preferred this approach rather than assigning arbitrary 
values. Also, the abundances for the strong-line elements did 
not change significantly. The best-fitting abundances are all 
reported in Table [T] 

To account for absorp tion, we used the Tubingen model 
(tbabs. [Wilms et al.ll2000T ). which should be regarded as an 
update of the common, and o lder, Wisconsin model (wabs, 
iMorrison fc McCamm"onlll983r i. A value of 4.0 • 10 cm 
was used for the Galactic absorption. 

The best-fitting values for the slope of the power-law 
component lie in the range 1.5-1.7, depending on the model 
(Table [T|, consistently with emission from X-ray binaries. 

The best-fitting \ 2 depends on the instrument whose 
data are analysed, on the plasma code (MEKAL or APEC), 
and on whether two lines possibly due to charge exchange 
(CE) emission (see Sect. 13 . 5f) are modelled. Not including 
the CE lines in the model, a switch from the Mekal- to the 
APEC-based model improves the fit: for the Mekal one, one 
gets, considering MOS+pn spectra, x 2 — 3172 with 2271 



at different temperatures) with a 6 order polynomial, and degrees of freedom (d.o.f.; x?cd = 1-40); for the APEC one, 
allows to vary the abundances. 

However, this model allows only temperatures with 
kT ^ 7 keV, causing the best-fitting DEM to have a sharp 
cut, clearly noticeable in Fig. [2] (left panel). Feeling that this 
limitation was somewhat artificial, we modified the available 
multi-temperature plasma models by enlarging the allocated 
energy range, so that the high-temperature part of the DEM 
could be better sampled (Fig. [2] right panel). The improve- 



ment in x was however not significant, indicating that the 



X = 2905 (Xrod = 1-28). A visual inspection of the residu- 
als confirms that the APEC model provides a slightly better 
representation of the data. On the other hand, after includ- 
ing the CE lines, the difference in x? g rea tly reduces, and 
both models yield similar values for the x? e d ■ The equivalent 
widths of the CE lines differ, between the Mekal and APEC 
models, by a factor of ~ 2-3. The spectra, with best-fitting 
models and x 2 values, are shown in Fig. [3j A comparison of 
the residuals is shown in Fig. [4] 
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Figure 2. The best-fitting DEM for the central regions of M82 by using the standard c6pvmkl model (left panel) and its enlarged energy 
array version (right panel). The total spectrum is computed by adding several components, one for each sampling point shown here. 



3.5 Charge exchange emission 

Two spectral lines which are not accounted for by the Mekal 
and APEC models, showed prominently in the residuals 
around 16 A (0.78 keV) and 10 A (1.2 keV); they may be 
best seen in the uppermost panel in Fig. [4] The best-fitting 
wavelengths and equivalent widts derived by parametrizing 
the lines with Gaussian models are reported in Table [T] 

The significance for all these lines when fitting MOS, 
pn, and joint MOS+pn data can be put at a confidence 
level > 99.97%. To esti mate the confidence level, we used 
Montecarlo simulations IjProtassov et al.ll2002l ) which can be 
briefly outlined as follows: 

(i) the model is fit to the data both with and without the 
line component ('line-model' and 'null-model', respectively); 

(ii) the variable 



Xnull ' 



Xlinc 



(1) 



(|Bevington fe Robinsonlll992h is calculated for the data; 

(iii) N sets of fake data files and fake backgrounds are 
generated according to the null-model. 

(iv) the fake data are fit with both the null-model and 
the line-model, and the F statistic as defined in Eq. ([TJ is 
calculated for each of them; 

(v) if M is the number of fake data sets for which F^ke > 
-Fdata, the confidence level is defined as 1 — M/N. 

These lines may be considered real also because: i) the 
effective areas in all EPIC cameras in the considered energy 
range are smooth functions of the wavelength and do not 
show any feature around the considered wavelengths; and 
ii) no background feature is present near these wavelengths. 
Moreover, the count rate of the central region of M82 is 
about 100 times larger than the background. 

We first looked for identifications in the Chianti 
database of lin e emission from astrophysical plasma 
l|Dere et alj|200ll ). The best candidates for the ~ 10 A line 
would be an Fe XVII line at 10.1340 A or a Ne X line at 
10.24 A, while for the ~ 16 A line the best candidate might 



be an Fe XVIII line at 15.83 A. However, the main prob- 
lem with this interpretation is that it is not clear why these 
lines should not be accounted for by the Mekal and APEC 
models. 

A more plausible hypothesis is that these lines, or at 
least the ~ 10 A one, are due to charge exchange (CE) 
emission, occurring at the interface between the hot wind 
and clouds of cold neutral gas within the galaxy itself. In 
the CE framework, ions from the wind diffuse through the 
interface into the cold gas, where electrons are transferred 
from the neutrals to the ions. Both the resulting ions can be 
highly excited and consequently re-arrange their electrons 
by emitting photons in the extreme ultraviolet and X-ray 
part of the spectrum. Only line emission occurs (no contin- 
uum), and the photon emission rate is proportional to the 
wind ion flux. CE emission has long been know n to occur 
in (optical) nebular spectra l)Chamberlainlll956h and it has 
been proposed to account for th e X-ray emission from the 
comet Hyak utake ([C ravens 1997), the solar system planets 
(e.g., Mars: iDennerl" et al. 2006), and it might also play a 
non-negligible role in explaining the emission from g alactic 
winds interacting with dense clouds l|Lallementll2004( ). 

In such a CE scenario, the ~ 10 A line may be iden- 
tified with emission from neutral Mg at 9.92 A (1.25 keV), 
and this would be consistent with t he presence of olivin e and 
pyroxene grams m the cold clouds l|Diuric et alJ l2005). The 
~ 16 A line, on the other hand, might be emission from ex- 
cited O VIII ions also resulting from CE interactions, where 
the electron is decaying from the n = 3 to the n = 1 level 
(but see also Sect 14. 2[) . 

iDiuric et al.l (|2005l ) suggest that, given the grain com- 
position cited above, the Mg line should be companied by a 
Si line at 1.73 keV. The addition of this line to the MOS+pn 
spectra is significant, at the 99.97% level, for the APEC 
model only, and does not remarkably change the Xred- It 
is, however, not significant for the MEKAL model (98.9% 
level). The equivalent width would be ~ 3.9 eV (APEC) or 
~ 4.6 eV (MEKAL). These values can be compared to the 
equivalent widths for the other two lines, which are found 
in the ranges 7-24 eV (see Table [1] for the details). 
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Figure 3. EPIC spectra of the central (2') region of M82. Upper panel: spectra with models. Lower panel: residuals, plotted in terms 
of sigmas with error bars of size one. A: observed spectrum (crosses) and model (solid line through the data points) convolved with the 
instrumental response of the pn. B: same, for the MOS1 (black) and MOS2 (red); these two cameras give nearly identical results so that 
the curves and points overlap. C and E: Gaussian line components of the model, for the pn. D and F: same, for the MOS1 and MOS2. 
G: absorbed power-law component of the model, for the pn. H: same, for the MOS1 and MOS2. I: thermal plasma (MEKAL) component 
of the model, for the pn. J: same, for the MOS1 and MOS2. Similar results are obtained with the APEC model (see also Fig. [4}. 



4 LARGE FIELD OF VIEW, HIGH 

RESOLUTION SPECTROSCOPY: RGS DATA 

The RGS instrument l|den Herder et all l200ll ) is a slitless 
grating spectrometer which allows, thanks to its better en- 
ergy resolution (E/ AE ~ 300, while the EPIC pn camera 
has E/AE ~ 7 in the overlapping energy range), a more 
precise identification of the emission lines and of their prop- 
erties, including their individual redshifts and — at least for 
broad lines — profiles. However, it is sensitive also to ex- 
tended emission on scales larger than ~ 1', and M82 was 
observed with the outflow extent along the RGS dispersion 
axis, so that for our purposes the line profiles are dominated 
by instrumental broadening. Thus we shall not discuss the 
line profiles in the following. 

The RGS spectra for the inner region were extracted 
with the rgsproc tool, considering only events within 90% 
of the PSF width (while the PSF width is actually energy- 
dependent, this corresponds roughly to ~ 1') in the cross- 
dispersion direction. On the dispersion direction, the in- 
coming photons are only limited by the instrument vi- 
gnetting. The alignment of the RGS cross-dispersion direc- 
tion is shown in Fig. [T] as the blue dashed lines. The back- 
ground spectra were taken from the blank sky observations 
using the rgsbkgmodel tool. Following the current calibra- 



tion documentation^]], systematic errors for an amount of 5% 
in the 9-12 A range, and of 8% in the 24.5-31 A range, have 
been allowed in the fitting procedure. 

The biggest problem in reducing the RGS data is the 
resolution degradation due to source extent. The rgsxsrc 
convolution model contained in the XSPEC package may be 
used to account for the source extent: it uses an image of the 
source to derive the line broadening function, assuming that 
the image profile is the same at all energies. Unfortunately, 
this is not true for M82 mainly because of the heavy ab- 
sorption to which its central regions are subject. To mitigate 
this problem, we used separate input images for the 6-18 A 
and 18-30 A bands, and for narrower bands when appro- 
priate. A source of uncertainty is the aperture size used by 
rgsxsrc to derive the broadening function. While a small 
aperture may introduce artefacts in the spectrum due to 
parts of the source not accounted by rgsxsrc yet contribut- 
ing to the spectrum, a large aperture may include too much 
background thus biasing the source profile. We found that 
an aperture of ~ 4' gives the best trade-off. 

In the fitting procedure used for the RGS data, the point 
sources contribution to the spectra is taken into account 



1 The CAL-TN-0030-3-0 document, available on the XMM web 
pages. 
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M. mospn 


M. mos 


M. pn 


M. mospn no-lines 


A. mospn 


A. mospn no-lines 


A^.th.MOS (10 22 cm" 2 ) 
JVff.th.PN (10 22 cm- 2 ) 


0.45 ±0.01 
0.44 ±0.01 


0.47 ±0.01 


0.44 ±0.01 


0.45 ± 0.01 
0.45 ± 0.01 


0.42 ±0.01 
0.41 ±0.01 


0.41 ± 0.01 
0.40 ±0.01 


O 

Ne[,Na] 
Mg,Al 
Si 
S 

Ar,Ca 

Fe,Ni 

norm th 


0.28 ±0.02 
0.42 ± 0.02 
1.00 ±0.02 
1.05 ±0.03 
1.05 ±0.05 
1.57 ±0.30 
0.33 ±0.01 
2.3 ■ 10" 4 


0.19 ±0.04 
0.49 ±0.03 
1.01 ±0.03 
1.16 ±0.04 
1.29 ±0.08 
1.83 ±0.47 
0.31 ±0.01 
2.4 ■ 10" 4 


0.29 ± 0.03 
0.35 ± 0.03 
0.83 ± 0.02 
1.01 ± 0.04 
0.98 ± 0.08 
0.85 ± 0.42 
0.33 ± 0.01 
2.4 • 10~ 4 


0.19 ± 0.02 
0.27 ±0.02 
0.72 ± 0.02 
0.93 ± 0.02 
0.98 ± 0.05 
0.74 ±0.25 
0.28 ± 0.01 
2.8 ■ lO" 4 


0.15 ±0.02 
0.46 ± 0.02 
0.77 ±0.02 
1.03 ±0.03 
1.18 ±0.06 
1.07 ±0.28 
0.30 ±0.01 
2.7- 10~ 4 


0.13 ± 0.01 
0.38 ± 0.02 
0.70 ±0.02 
0.96 ± 0.02 
1.10 ±0.06 
0.48 ± 0.26 
0.29 ± 0.01 
2.6 • 10~ 4 


JVff.pow (10 22 cm" 2 ) 

r 

norm pow 


2.26 ±0.09 
1.51 ±0.03 
1.2 ■ 10~ 3 


5.4 ±0.5 
1.48 ±0.09 
8.3 ■ 10~ 4 


3.5 ±0.4 
1.7 ±0.1 
1.2 • 10~ 3 


5.25 ± 0.27 
1.62 ± 0.05 
1.1 ■ 10~ 3 


4.70 ±0.25 
1.59 ±0.05 
1.0 • 10~ 3 


5.27 ±0.25 
1.63 ± 0.05 
1.2 • 10~ 3 


x 2 


2743 (2267) 
1.21 


1101 (881) 
1.25 


1471 (1365) 
1.08 


3172 (2271) 
1.40 


2779 (2267) 
1.23 


2905 (2271) 
1.28 


line energy (keV) 
line wavelength (A) 
line eq. w. (eV) 


1.234 ±0.011 
10.05 ±0.09 
22.6 


l.24ol°;?3 6 

15.8 


i 91Q +0.007 
1 - zla -0 .005 
10 17+ 004 

15.8 




i 999 + 0.004 

i.zzz_ Q 0Q2 
10 15+ 002 

iu.io_ 04 
7.74 




line energy (keV) 
line wavelength (A) 
line eq. w. (eV) 


„ +0.002 
"■III -0.006 

io.ao_ 03 
22.4 


n 7 on+0.003 
U. .007 

15.82+^ 4 
31.7 


n 771 +0.008 

u - ' ' 1 -o.oo6 
16.08±°;£ 
16.4 




0.767 ±0.006 
16.15 ±0.12 
13.6 





Table 1. Spectral fit parameters for the EPIC data in the inner region. All quoted errors are at 90% confidence level. The rows from top 
to bottom show: column density of absorbing material relative to the thermal component in MOS and PN data; chemical abundances 
of the thermal component relative to solar values; normalisation of the thermal component (units as in Eq|4]l; column density of the 
absorbing material relative to the power-law component; power-law photon index; normalisation of the power-law component (in units 
of photons s — 1 cm -2 kev - 1 at 1 keV); x 2 an d reduced x 2 ; equivalent widths of the CE lines. The quoted errors are referred to the 90% 
confidence level, and have been calculated with all other parameters in the spectral fit (DEM, abundances, Nh , power-law parameters) as 
frozen. If this assumption is relaxed, e.g. by thawing other element's abundances, the O, Ne and Fe errors relative to the MEKAL-mospn 
case would be ±0.03, ±0.02 and ±0.01 respectively. By relaxing it even more, e.g. by thawing also the six DEM parameters, one would 
get ±0.05, ±0.05 and ±0.03 respectively. 



by including the same power-law component discussed in 
Sect. [3] with fixed parameters. 



4.1 Thermal lines 

The spectra, with the best-fitting model and the residuals, 
are shown in Fig. [5] (6-18 A) and Fig. [7] (18-30 A), while 
the spectral results are listed in Table [5] The brightest lines 
between 6 and 10.5 A are from Si XIII (a blended triplet 
around 6.7 A), Mg XII (A8.42 A), and Mg XI (A9.17 A). 
Other elements with fainter emission in this region are Al, 
Na and Ne, while the 10.5-12.0 A region is dominated by 
L-shell Fe emission, mainly from Fe XII, Fe XIII and Fe 
XIV. The Ne X A12.1 A line is also prominent. At longer 
wavelengths many lines are present from Fe and Ni, plus 
some bright lines from Ne and O; the brightest ones are Fe 
XVII A15.0, 16.8-17.1 A; Fe XVIII A14.2, 14,4 A; Fe XIX 
A13.5 A), and from Ne IX (A13.4 A) and O VIII (A16.0 
A). Longwards of 18 A, the lines O VIII A19 A, O VII 
A21.6-22.1 A, N VII A24.8 A and a hint of C VI A28.5 A 
are detecte d. Line identifica tions and wavelengths are from 
CHIANTI (|Dere et al.ll200ll . and references therein). 

The O lines, because of the possible contribution from 
the CE mechanism and because a different best-fitting red- 
shift is needed for O VIII A 19 A than for the lines at shorter 
wavelengths, will be discussed in Sects. 14.21 Thus, for the 



dependence on energy of the line profiles and also because 
the N VII and C VI lines have much fewer counts than all 
other lines, it is best to restrict the interval to analyse to 
the 6-18 A region. This restriction does not affect most ele- 
ments (i.e. Ne, Mg, Al, Si, Fe and Ni) whose lines are weak, 
if present at all, in the 18-30 A region. 

The best-fitting values for the abundances provided by 
the two codes are in agreement within errors for all elements. 
The APEC model gives the best-fitting to the observed 6- 
18 A spectrum, with Xred ~ 1-15. It was not possible to get 
a good fit by using the MEKAL model (the best one has 
Xr C d ~ 1-38). The best-fitting abundances from both models 
are shown in Tabled and the APEC-derived ones are also 
displayed as the red horizontal lines in Fig. I12[) . 

The restricted wavelength interval 12-18 A looks like 
the best one in which to test our assumption that the Fe 
and Ni have equal abundances, both because of fit good- 
ness (x 2 ~ 1.05) and of the relative paucity of lines from 
other elements. By using an F-test, the improvement in x? G d 
resulting from thawing the Ni abundances was not signifi- 
cant (improvement only at 75% confidence level); this was 
achieved by keeping all other parameters except Fe frozen; 
thawing O and Ne further reduced the significance of the x 2 
improvement. 
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Figure 4. Comparison of residuals of EPIC spectra of the centre of M82 from different models. From top to bottom: Mekal without CE 
lines; APEC without CE lines; Mekal with CE lines; APEC with CE lines. The x 2 cd values are 1.40, 1.28, 1.21 and 1.23, respectively 
(see also Table[T]l. Black points: MOS1; red: MOS2; green: pn. 



The N VII and C VI lines 

Because of the strong absorption, few photon counts are 
available in the 24-30 A region where the N VII and C VI 
lines are present. No significant constraint may be put on 
the C abundance, which would mainly be obtained from the 
C VI line at 28.46 A where the absorption is maximum. 
An estimate of the N abundance may be placed at values 
larger than solar, albeit with rather large errors. The random 
errors are reported in Table [2] As to the systematic ones, 
one could make the following consideration. The best-fitting 
abundance of N in Table [2] was calculated by using an input 
image for rgsxsrc in the 24-26 A band, and the redshift 
as a free parameter. Had we included the 24-26 A region 
along the 6-18 A one, with the same 6-18 A input image for 
rgsxsrc, and assigned to the N VII line the same redshift of 
the other considered ones, then the best-fitting N abundance 
would have been 1.0 ± 0.2. 

4.2 The O lines 

Two lines from O VIII (A16 A and A19 A) and one triplet 
from O VII (21.6-22.1 A) are present in the RGS spectrum. 
It is not possible to fit all lines together, either with the 
APEC or with the MEKAL models, because these models 





Mekal 


Apec 


N H (10 22 cm~ 2 ) 


0.21 ±0.01 


0.15 ± 0.01 


N 




3.2 ± 1.6 


O (16A) 


1.31 ±0.16 


0.62 ± 0.10 


Ne[,Na] 


0.96 ±0.06 


0.95 ± 0.07 


Mg 


1.69 ±0.12 


1.68 ± 0.10 


Al 


1.9 ±1.3 


3.1 ± 1.2 


Si 


2.49 ±0.23 


2.64 ±0.23 


Fe,Ni 


0.51 ±0.02 


0.55 ± 0.02 


normth 


2.28 • 10~ 4 


2.15 ■ 10~ 4 


K 


1559.7 (1133) 


1303.4 (1133) 


•*4cd 


1.38 


1.15 


O (19A) 


0.87 ±0.03 


1.41 ± 0.05 




284.1 (148) 


181.3 (148) 




1.9 


1.22 


eq. w. of Mg line (eV) 


13.1 


12.4 


eq. w. of Si line (eV) 


43.7 


43.3 



Table 2. Chemical abundances from spectral fits to RGS data. 
The quoted errors are referred to the 90% confidence level, and 
have been calculated with all other parameters in the spectral fit 
(DEM, abundances, Nh, power-law parameters) as frozen. 
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Figure 5. RGS spectra (upper panel) and residuals (lower panel), 6-18 A. Black data points, with markers: RGS1; red data points, 
without markers: RGS2. The most intense lines are labelled with their identifications. In this and in the following pictures about RGS 
data, the points with very low counts looking like wells in the spectrum represent the gaps between two consecutive RGS CCDs or bad 
columns. 



APEC MEKAL 




channel wavelength (Angstrom) channel wavelength (Angstrom) 

Figure 6. Comparison of the APEC (left) and MEKAL (right) modellings of the O VIII A16 A line. The goodness-of-fit is x^ od = 1-0 
for the left panel and X^cd = ^ f° r the right one. A different modelling of the Fe XVIII lines at 15.4 A, 15.6 A and 15.8 A is probably 
the origin of the discrepancies. Crosses: data points. Solid lines: models. Black: RGS1. Red: RGS2. 



A deep X-ray observation of M82 9 



assume a common redshift for all lines, which is not the case 
for O. Thus these lines will be analysed separately. 

An additional complication may be the possible pres- 
ence of CE emission, that should company the Mg and Si 
emission if the ions impacting on the dust contain O. The 
relative strengths of the three components of the O triplet 
arising from CE have recently be en measured for the firs t 
time in laboratory experiments ( Beiersdorfer et all [2003), 
finding that the dominant contribution comes from the for- 
bidden line. This effect has also been observed astronom- 
ically in J upiter's auroras, albeit with large errors on the 
line fluxes (|Branduardi-Ravmont et al.l 2006). Experiments 
on CE emission from the Fe XXV triplet have also con- 
firmed the predominance of the intercombi nation and for- 
bidde n lines relative to the resonance one (|Wargelin et al.l 
l2005h . 



4.2.1 The O VIII A 16 A line 

The O VIII A 16 A line can be fit together the other lines 
m the 6-18 A range; the abundances derived from this line 
are shown in Table [2] along with the results from the other 
elements. A juxtaposition of close-up images of the O VIII 
A16 A with both the APEC and MEKAL models is shown in 
Fig. [6] The goodness's of fit in the 15.75-16.45 A interval are 
X% d = 1.0 and x? cd = 1-6 for the APEC and MEKAL mod- 
els, respectively. From the comparison, one may attribute 
the deviations to a different modelling of Fe XVIII lines at 
15.4 A, 15.6 A and 15.8 A. The APEC-derived value for the 
O abundance seems therefore more robust. 



The O VIII A19 A line 

The O VIII line A19 A formally needs a different redshift 
than the lines at shorter wavelength to be fitted. The best- 
fitting redshifts are 9.6 • 10" 4 and for the APEC and 
MEKAL models, respectively. However, the relative shifts 
between O VIII A19 A and the lines at shorter wavelengths 
are of the same magnitude of the systematic uncertainties 
(see Sect. 14.4)) . The uncertainties are mainly due to the de- 
pendency of the line broadening function on the energy. The 
amount of absorption makes O lines especially sensitive to 
this effect, so that using a narrow band image as the input 
for rgsxsrc gives the best results. 

The O abundances extracted from this line, with all 
other parameters kept the same as for the 6-18 A spectrum, 
are shown in Table [5] While the APEC model gives an ac- 
ceptable fit (Xicd ~ 1-2) by using the same parameters for 
the 6-18 A spectrum, the MEKAL model needs a differ- 
ent value for the column density or for the normalisation to 
achieve a good fit. For the MEKAL model, had we thawed 
the column density, we would have got Nu ~ (4.5±0.1) • 10 21 
cm -2 , O/O ~ 3.6±0.1 and xLj ~ 1-2. Had we thawed the 
normalisation, a similarly high abundance of O would have 
been obtained. 



The O VII triplet 

In the 20.5-23.0 A the only strong lines are the He-like O 
VII triplet, at 21.6 A, 21.8 A and 22.1 A, known as the res- 
onance, intercombination and forbidden lines, respectively. 



We shall retain this nomenclature hereafter. One difficulty 
in fitting the O VII triplet is that all the lines are blended. 

The RGS1 spectrum in the O VII region (RGS2 is not 
working at these wavelengths) is reproduced neither by the 
APEC nor by the MEKAL model. Thus in the following we 
will try and fit the spectrum with Gaussian lines. The fol- 
lowing line ratios are used in the literature, which have been 
show n to be sensitive to the electron density and tempera- 
ture (|Gabriel fe Jordanll 19721 ): 



%) = f/i 
G(T e ) = (/ + i)/r 



(2) 
(3) 



where /, i, and r represent the fluxes in the forbidden, in- 
tercombin ation and resonant line re spectively. 

From iBeiersdorfer et ail (|2003f) one would get, for CE 



emissi on, R ~ 3 and G ~ 2.2; iBranduardi-Ravmont et al.l 
J2006I) gives values which are consistent, within errors, with 
IBeiersdorfer et all 's results. 

For a thermal plasma dominated by collisional ionisa- 
tion (CI), the ranges for these parameters are 2 ■ 10~ 2 < 
R < 4 (considering densities ranging from 10 12 cm -3 to 
1 cm -3 , respectively) and 0.6 < G < 1.2 (temperatures 
from ~ 0.35 keV to ~ 0.1 keV, respectiveljQ) ; however the 
presence of photo-ionisation may lower R mimicking a high 
density plasma; for a hybrid (coll isional+photo) plasm a one 
should consider 2- 10~ 4 < R < 4 jPorquet et alj|200ll ). Thus 
it may be difficult to distinguish between a high density, low 
temperature collisional pl asma and CE emissio n. In the fol- 
lowing, we always refer to lPorquet et al.l l|200ll ) for the links 
among R, G and the physical parameters of the plasma. 

It is possible to obtain an acceptable (x? c d ~ 1-2) fit to 
the spectrum by using a single triplet of lines, provided that 
some additional broadening is present beyond the ones ac- 
counted for by the instrumental response and rgsxsrc; this 
may be parametrized as a dispersion a ~ 3.7 ± 1.3 eV for 
each Gaussian component of the triplet. If this broadening 
were of thermal origin, the corresponding velocity would be 
~ 8 km s _1 . The line intensity ratios are marginally consis- 
tent with CE but not consistent with plasma emission, as 
may be seen in the confidence contours plotted in Fig. [8] 
It is not clear why only O VII needs additional broadening, 
with respect to all other lines; while it may be that multiple 
components (for instance, CE and plasma) are present at 
different velocities, the data resolution is not high enough 
to disentangle this issue. 



4.3 Lines with possible contribution from 
Charge-Exchange: Mg, Si 

Some residuals around 7.2 A could be identified with a CE 
line from neutral Si, with a best-fitting energy within errors 
from what expected (rest fra me wavelength: 7.1280 A, or 
1.7394 keV. iDiuric et al.ll2005l ) . The significance of this line 
can be estimated, with the Montecarlo simulations described 
in Sect. l3.5l and using the APEC model, at a confidence level 
> 99.99%. 

Another possible identification for this line might be Al 

2 While temperatures warmer than 0.35 keV may also be plausi- 
ble, we did not consider this possibility in the fit, since for higher 
T c the flux in the whole O VII triplet decreases. 
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channel wavelength (Al channel wavelength (Al channel wavelength (A) 

Figure 7. RGS spectra and residuals in the 18—30 A region. From left to right: spectrum and residuals of the O VIII A19 A line, fitted 
with the APEC model; same, for the O VII triplet with best-fitting model; same, for the N VII and C VI lines fitted with the APEC 
model. 
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Figure 8. Confidence contours (68.3%, 90% and 99%) for the R 
and G parameters (which describe the line intensity ratios) for 
the O VII triplet. The spectrum is marginally consistent with CE 
emission (shown as the grey circle; the size of the circle is not 
representative of the uncertainties), but is not consistent with 
plasma emission (hashed rectangle). 



XIII A7.17 A, which however is not allowed by the APEC 
model because an enhancement of this line due to a larger 
Al abundance should be companied by a similar increase in 
the Al XII triplet at 7.76-7.87 A, which is not observed. In 
order to reproduce the Al XIII/ Al XII intensity ratio, the 
plasma temperature should be kT ~ 2 keV, which is at odds 
with the DEM structure derived from the other elements. 

The Mg CE line (referred to as 'the ~ 10 A' line 
in Sect. 13. 5[) is detected only at a lower confidence level 
(99.8%). 



4.4 Comparison of best-fitting velocities 

In Fig. [9] we show a comparison of the different ranges for 
the best-fitting velocities of the X-ray plasma in the 6-18 
A range, the O VIII A19 A line, the O VII A21.6-22.1 A 
triplet, and the N VII A25 A line. The lines in the 6-18 A 
interval formally have a larger velocity than the others at 
longer wavelengths. 

In the following we discuss possible systematic uncer- 
tainties which may be responsible of the observed discrep- 
ancy between 'short-' and 'long'-wavelength lines. 

• We recalculated the Xicd distribution for the redshift of 
the N VII line using a 6-18 A image, in order to maximize 
the distance between the wavelengths of the line and of the 
reference image. The result is shown in grey in Fig. [9] with 
the label 'N VII check rgsxsrc'. We found that in this case 
the redshift of the N VII line becomes slightly larger, but 
it is still not consistent with the redshift of the lines in the 
6-18 A range. 

• Using a 6-12 A image as input for rgsxsrc while study- 
ing the O VIII A19 A line results in a much worse fit 
(Xred ~ 1-4) than using an image in the proper band 
(Xred ~ 1-0). Allowing additional broadening does not make 
the fit better. 

• At variance, splitting the 6-18 A range in two smaller 
intervals (6-12 A and 12-18 A) did not show any difference 
in the best-fitting velocities, with both cases confirming the 
6-18 A value. 

The shift between the 'long-' and 'short'-wavelength 
lines amounts to Az ~ 3- 10 -3 ; considering an average wave- 
length of 15 A, this corresponds to a AA = A Az ~ 45 mA, 
hence to a spatial separation of ~ 20", if due to an angular 
offset of the two emitting regions within the RGS aperture. 
This distance is shorter than the aperture considered for 
rgsxsrc, and is comparable to the size of the region where 
there is the maximum foreground absorption. Thus the ve- 
locity shift may well be possibly attributed to the apparent 
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Figure 9. Comparison of best-fitting velocities, relative to the bulk motion of M82 (z = 6.77 ■ 10~ 4 . lde Vaueouleurs et al,lll99lh . for 
different lines and wavelength intervals in the RGS spectra. The crosses show the best-fitting values, while the boxes represent the 68.3% 
confidence interval, and the solid- and dotted-line error bars show the 90% and 99% confidence intervals respectively. The grey whisker 
labelled 'check rgsxsrc' shows the maximum possible systematic error arising from an incorrect image supplied to rgsxsrc, estimated for 
the N VII A25 A line using a 6— 12 A image. 



spatial separation of the regions emitting at long and short 
wavelengths. 



5 SPECTROSCOPY OF THE OUTFLOW 

To study the different properties of the hot gas as it flows 
and is heated from the central starburst towards the inter- 
galactic space, we divided both outflows in five regions each. 
Each region has a rectangular shape, the larger side being 
parallel to the galaxy major axis. The regions are numbered 
from Nl to N5, and from SI to S5 with increasing height 
above the galactic plane as shown in Fig. [T] in the figures 
we adopt the convention of assigning negative heights to the 
northern regions. The regions have different sizes, chosen 
mainly to guarantee enough photon counts and a sufficiently 
good signal to noise ratio to perform a meaningful spectral 
analysis in each area. The spectra were extracted from the 
MOS and pn data, and fitted with the same models used for 
the central region, the only difference being point sources, 
which in the outflow are sufficiently apart to be excluded 
when extracting spectra. The positions of point sources were 
determined by running the PWXDETECT wavelet-based 
detection code on the MOS1, MOS2 and pn images (the 
algorithm is described in lDamiani et al" 1997). An atlas of 
the spectra with best-fitting models is shown in Fig. 1101 The 
APEC model was used; no significant difference was found 
using the Mekal one. 

The same multi-temperature structure is found in every 
region. By fitting the DEM polynomial with two Gaussian 
components, we found that the centre of the colder com- 
ponent moves slightly from 0.53 keV in the centre of M82 
to 0.33 and 0.37 in the outer regions (N5 and S5, respec- 
tively) . The hotter component does not significantly change 



its temperature. Fig. [TT] shows the centres and FWHM of 
the best-fitting Gaussians. 

The best-fitting chemical abundances from EPIC and 
RGS data are shown in Fig. [12] along with results from 
infra red observations for the central regions l|Origlia et al.l 
2004). Only one region (SI) required in its spectrum some 
absorption in excess of the Galactic value, probably due to 
the M82 disc partially coveri ng the southern ou tflow. We 
confirm our previous finding (|Origlia et alj 12004). that the 
lighter Q-elements are more concentrated in the outflow than 
in the centre. This effect is larger for elements with lower 
atomic mass, becomes less evident for Si and reverses for S. 
The centre/outskirt abundance ratio in the centre is about 
~ 1/10 for O and Ne. Fe is also more concentrated in the 
outflow. 

The abundances ratios, shown in Fig. 1131 have smaller 
variations, and present different trends for light and massive 
elements: while the O/Fe and Ne/Fe ratios are lower in the 
centre than in the outskirts, the opposite holds for Si/Fe 
and S/Fe, with Mg/Fe being an intermediate case showing 
no variation. The scatter between values for centre and out- 
skirts is a factor of < 3. 

The physical parameters of the plasma may be obtained 
from the temperature (T) and normalisation (A) of the 
model, with some assumptions about the volume (V) and 
filling factor (/). The normalisation of plasma models in 
XSPEC is defined as 



47r(D A (l + z))~ 2 ■ J n c n n dV 

using c.g.s. units and where Da is the angular distance, and 
n e and nu are the electron and hydrogen volume densities. 
One usually assumes n c ~ nu, and defines the Emission 
Integral (EI) as EI = J n 2 dV. From this, one can obtain 
the density n c ~ (EI/Vf) 1/2 , the pressure p ~ 2n c kT, 
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Figure 10. EPIC spectra for increasing height on the galactic plane. Green points (upper data series): pn. Black and red points (lower 
data series): MOS1 and MOS2, respectively. 
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Figure 12. Variation of chemical abundances with increasing height on the galactic plane. Black: abundances from X-ray MOS and pn 
data. Blue: abundances from X-ray RGS data (due to the characteristics of the RGS spectrometer, they represent space-averaged values). 
Red: abundances from infrared data (corresponding to red supergiant stars in the galaxy central region). Negative values of distance 
refer to the north direction, positive values to south. 



the gaseous mass M ~ n e m p Vf, the thermal energy con- 
tent E ~ 3n e kTV f , and the radiative cooling time t ~ 
3fcT/(An c ), where m p is the proton mass, and A is the emis- 
sivity of the gas. Since A depends on th e temperature, its val- 
ues w ere interpolated from Table 6 of ISutherland fc Dopital 
l|l993t l (A is of the order of - 3 • 10" 23 erg cm 3 s" 1 , valid 
for solar metallicity and the temperatures of interest). 

The parameters were calculated using the best-fitting 
normalisations for each region. The used temperatures are 
the best-fittings to the low-energy components of the DEMs 
(see lower panel in Fig. For each of the rectangular N 



and S regions (Fig.[TJ, the volume was assumed to be that of 
a cylinder with base diameter and height equal to the major 
and minor sides of the rectangle, respectively. Normalisa- 
tions, volumes and physical parameters of the plasma are 
shown in Table [3] The values for the galaxy centre are also 
shown in the same table, assuming a cylindrical volume with 
base radius 32" and height 120". From the Table, one may 
see that the gas density and pressure decrease by a factor 
of ~ 10 from the centre to the outskirts, while the cooling 
time increases. 
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N5 


N4 


N3 


N2 


Nl 


centre 


SI 


S2 


S3 


S4 


S5 


Normalisation (10 5 ) 





14 


0.68 


1.6 


1.6 


4.0 


27 


5.3 


2.4 


0.77 


0.14 


0.18 


Volume (kpc 3 ) 





61 


1.2 


1.1 


0.43 


0.41 


1.4 


0.24 


0.35 


0.26 


0.21 


0.22 


Density (10 -3 cm -3 ) 


3.4/ 


i 


5.4 


8.7 


14 


23 


32 


35 


19 


13 


6.0 


6.5 


Pressure (10 — 12 dine cm -2 ) 


3.6/ 


i 


6.1 


10 


19 


35 


54 


58 


31 


19 


8.3 


7.7 


Mass (10 5 Mq) 


0.52/ 




1.7 


2.4 


1.5 


2.3 


11 


2.0 


1.7 


0.81 


0.31 


0.37 


Energy (10 53 erg) 


0.99/ 




3.3 


5.2 


3.6 


6.3 


33 


6.0 


4.8 


2.2 


0.76 


0.77 


Energy density (10 — 12 erg cm -3 ) 


5.5/ 


H 


9.2 


15 


28 


52 


81 


87 


47 


28 


12 


12 


Cooling time (Myr) 


450/ 




310 


210 


160 


120 


100 


89 


160 


220 


390 


280 



Table 3. Physical parameters of the plasma across the different regions of the outflow. The dependencies on the filling factor / have 
been explicited for simplicity only in the first column, but they apply to all columns. 
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Figure 11. Variation of temperatures of the plasma along the 
M82 outflow. The lower and upper panels show the first (~ 0.5 
keV) and second (~ 7 keV) peak of the DEM, respectively. The 
error bars represent the FWHM of Gaussian curves fitted to the 
peaks of the polynomial distribution of the DEM. For instance, 
the two error bars at kpc and labelled 'ctr' (i.e. centre) are just 
a different representation of the right panel of Fig. [2] Negative 
values of distance refer to the north direction, positive values to 
south. 
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6 DISCUSSION 

6.1 Multi-temperature nature of the gas 

To find why the DEM has a double-peaked shape (Figs. [2] 
and we first checked if the hot component (the peak 
at ~ 5 keV) was really required to achieve a good fit. 
By truncating the DEM and leaving only the low energy 
(~ 0.5 keV) peak, no similar or better fit than those shown 
in Table [T] could be obtained, {\ 2 ~ 3028 with 2273 de- 
grees of freedom, Xrcd ~ 1-33, CE lines included). Then we 
checked if the hot component could have a non-thermal ori- 
gin. Since star forming galaxies are known as radio sources, 
bremsstrahlung emission should in principle arise from the 
same population of electrons from which the synchrotron 



Figure 13. Abundance ratios (X/Fe) observed in EPIC spectra of 
the outflow. Negative values of height refer to the north direction, 
positive values to south. 

emission is originated. Thus we developed a non-thermal 
bremsst rahlung model for XSPEC following the prescrip- 
tions in iBlumenthal fc Gouldl l|l970l) . By including such a 
component in the model, no improvement of the fit statis- 
tics was observed with respect to the previous case (i.e., no 
hot component at all). 

It is unlikely that the hot component is due to unre- 
solved point sources. By analysing the Chandra 3-9 keV 
spectra of the diffuse emission of M82, and by considering 
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the Lor N-Log S of point so urces in star forming galaxy, 
IStrickland fc Heckmanl (|2007l ) have shown that significant 
flux remains in that band after accounting for unresolved 
sources. 

Thus, assuming that both components were thermal, 
one might naively identify the two temperatures with plasma 
coming from different regions, but this idea would have a 
hard time explaining why the temperatures are about the 
same across all the outflow, and also why a similar pat- 
tern with same temperatures seems to be exhibited by other 
galaxies, e.g. NGC253. 

The temperature pattern exhibited by M82 is strikingly 
similar to that of the Galactic ridge (GR), which is an en- 
han cement of the X-ray background along the Galactic equa- 
tor (|Cooke et al.ll 19691 ), and whose explanations range from 
truly diffuse emission to the sum of many faint sources. Un- 
like the GR case, an interpretation based on the hypothesis 
of unresolved point sources does not seem viable for M82. 
The main evidence against is the morphology of the outflow, 
which is not associated with any stellar structure, such as 
the disc or the halo. 

The temperature pattern is not the only similarity be- 
tween M82 and the GR: another, important one is the en- 
ergy of the Fe K lines. By using a simple model (power-law 
plus Gaussian line) for the 4-8 keV EPfC spectrum, an Fe 
thermal line was detected! at 6.66 ± 0.02 keV (cr ^ 0.08, 
X? od ~ 622/568 ~ 1.09) with equivalent width ~ 100 eV 
(th us confirming the v alue originally found with BeppoSAX 
bv ICappi et al.lll999h . A similar value is found in Suzaku 
data (P. Ranalli, in prep.). 

For the GR, iKovama etlll (jl986T ) found a line energy 
of 6.71 ± 0.04 keV in observations with the Tenma satellite. 
iKaneda et all (|l997f ) reported a line energy of ~ 6.6 keV. 
The best-fitting line width also seem to agree (GR: a ~ 70 
eV; M82: a < 80eV). One should also note, however, that 
ITanakal ([2002) reported a higher energy for the Fe line in 
the GR (6.76 keV), and that a line centroid of ~ 6.66 keV 
is also consistent wit h a CE origin for the Fe XXV triplet 
l|Wargelin et alJl2005h . 

Anyway, the presence of hot gas with kT ~ 7 keV should 
put a significant number of Fe atoms in the highest ioni- 
sation states, namely He- and H-like, whose line emission 
happens at higher energies (~ 6.8 keV). Since the presence 
of such a line in M82 can be safely excluded, one might ask 
if the hot plasma component is indeed real or is just an 

artefact. 

An im portant insight has been provided bv lDogiel et al.l 

l|2002l ) and lMasai et al.l (|2002h . The authors argue that the 
bremsstrahlung emission of electrons from the regions in 
which particle acceleration occurs has a peculiar spectrum 
which mimics the double-temperature pattern found both 
in the GR and in M82. By considering in situ acceleration 
of thermalised electrons, they derive a broad transition re- 
gion between the thermal and non-thermal energy regions. 
The total X-ray spectrum resulting from both the thermal 
and transition regions of the electron energy spectrum may 
then be described, with good approximation, by the sum 
of two thermal plasmas whose temperatures are in the ra- 



3 This line is fully reproduced by the MEKAL and APEC models 
of Sect. |3] 



tio Thotter/Tcoider ~4-5. In this model, the actual temper- 
ature of the thermalised gas is the colder of the two (the 
authors propose a value of 0.3 keV for the GR). This would 
then naturally explain the measured energy of the Fe K 
line. Another consequence of this model would be the pres- 
ence of synchrotro n emission from the outflow , which in M82 
is already known l|Seaauist fc Qdegard|[l99lh . However, the 
Tkotter/Tcoidor ratios observed in M82 are larger than the 
predicted one by a factor 3-4 (see Fig. 1 1 1 p . While the model 
looks promising, some more work is still needed to check if 
the discrepancy in the temperature ratio can be reconciled. 

6.2 Vertical abundance gradients 

In a somewhat naive scheme, a cavity of hot enriched gas 
is created in the ISM when the first supernova explosions 
begin. The hot gas which happens to be at the boundary 
between the cavity and the cold ISM may either cool and 
mix, or stay hot and unmixed. Let us suppose it does not mix 
significantly until the hot gas bubble breaks out of the galaxy 
plane and the outflow is formed. Since the more massive a 
SN progenitors is, the more metals are produced, the outer 
regions of the outflow should have larger abundances than 
the inner one. 

To set this idea in a quantitative framework, the rela- 
tionships between the lifetime of the SN progenitor and the 
mass of material releas ed after the explosion are shown in 
Fig. 1141 with data from IWooslev fc Weaver! (1 19951. h e reafte r 
WW95) and evolutionary tracks from Bressan et al.1 l|l993f ). 
Among the several models present in WW95 we have con- 
sidered a solar value for metallicity, since this meets most 
closely the IR and X-ray values. In Fig. [15] we also show the 
X/Fe mass ratio from WW95. 

To locate the proper time-scale, it is important to con- 
sider the age of the starburst in M82. iMatsushita et all 
(2000) report from CO measurements the presence of a 
molecular super-bubble of size ~ 120 x 240 pc in the cen- 
tral regions of M82, expanding at ~ 100 km s _1 , which lead 
them to derive an age of 1-2 Myr. However, this bubble 
seems to have not yet fully broken, so it may represent a 
starburst episode succe ssive to the one(s) that p roduced the 
X-ray emitting outflow. [Stricklan d et al.1 (|l997l ) report from 
ROSAT observations of the X-ray halo that the age may be 
10 million years. Also. lSilva et al.1 (|200lh investigated the ra- 
dio to X-ray spectral energy distribution in M82 and found 
a value for the starb u rst ag e of < 25 Myr. The models in 
IStrickland fc Stevens! (|2000T ) suggest that the time required 
for the starburst in M82 to break out of the ISM is about 
5 million years. Finall y, the presence of r ed supergiant stars 
in the centre of M82 ( Origlia et al.ll2o"o"4h puts a lower limit 
to the starburst age at ~ 7 Myr. 

Looking at Fig. [14] one may find a global decrease of 
the production of O, Ne, Mg, Si, S for increasing progeni- 
tor lifetime (or, equivalently, for decreasing mass). In other 
words, the first stars to explode are also the ones with larger 
yields. The Fe production, on the other hand, has a less pro- 
nounced decrease, with a secondary maximum at ~ 18 Myr; 
it also has large uncertainties for short-lived stars. It may be 
tempting to look again at Fig. 1121 under the hypothesis that 
the emission further away from the galaxy centre traces the 
ejecta from more massive and short-lived stars. This picture 
may be checked also by comparing the observed and the- 
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Mass (M e ) 

Figure 14. Masses (in Mq) of ejected material from type-II SN 
explosions of massive stars with solar metalli city, as functions of 
the li fetime in the main sequence (data from IWooslev fc Weaver! 
1 19951) . The (black) solid line refers to calculations made using the 
standard value of l.e ■ 10 51 erg for the final kinetic energy of the 
ejecta (model A in WW95). The (red) dashed and (blue) dot- 
dashed lines refer to more energetic explosions (models B and C 
in WW95) which affect only the mo st massive progenito rs. For 
clarity, a mass scale is also displayed feressan et al.l l 1993). 



oretical abundance ratios (Figs. 1131 and [T5 |) . However, sev- 
eral pitfalls may hamper this interpretation. For instance, 
we have ignored here any role of the IMF or of the mixing; 
this is only true insofar as one assumes that the outflow gas 
is produced in an instantaneous starburst, and then moves 
in bulk with a constant velocity and without running into 
other material, such as cold neutral gas in the disc or hot 
halo gas. To cope with this difficulties would require more 
assumptions about the history of the starburst-produced gas 
and of its outflow than it is advisable to discuss here. For in- 
stance, if the starburst is not instantaneous but rather has a 
finite duration, then one needs to consider the IMF, whose 
effect is to introduce a blurring in the yields distribution 
of Fig. 1141 Also, the mixing of the starburst-produced gas 
with pre-existent halo gas might veer the observed chemical 
abundances ratios from resembling type-II SN yields towards 
type-la SN yields. 

Another possibility is that the M82 wind may be signif- 
icantly more mass-loaded in the galaxy inner regions than 
in the outer ones. It is known that the X-ray emission from 
the wind regions is not due to a volume-filling wind fluid, 
but rather it arises at the inte racting boundaries of the wind 
and of the ambient material (|Marcolini et alj|2005l . and ref- 
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Figure 15. Abundance ratios (Mx/My) of material ejected by 
type-II SN explosions of massive stars with solar metallicity, 
as functions of the lifet ime in the main sequence (data from 
IWooslev fc Weaver! Il995h . Colours, line styles and scales as in 
Fig.QH 



erences within). Thus the metallicity of the X-ray emitting 
gas results from the mixing of the wind with the ambient 
material. 

The wind/cloud interaction is related to the mass load- 
ing of the wind, i.e. whether the clouds survive to the wind 
encounter, or on the contrary if they are disrupted, their 
material is accelerated and becomes part of the wind. The 
lifetime of the clouds is not very well known, since it depends 
on how thermal conduction i s modelled at the wind/cloud 
border. iMarcolini etafl (|2005T ) estimate for a cloud of radius 
R = 15 pc a lifetime in the 0.9-12.8 Myr interval, depend- 
ing on the conduction model. It is also known, however, that 
it scales with R 2 , R being the cloud radius, so that larger 
clouds may have much longer lifetimes. 

Thus, recalling that the starburst age in M82 is < 25 
Myr, it is possible that in the inner regions on M82 there 
are still some clouds not yet disrupted. If the total surface of 
the clouds which is exposed to the wind has increased with 
time, e.g. by fragmenting, than it may be possible that mass 
loading is becoming more and more important. Provided 
that the supernova-originated wind is enriched with respect 
to the ambient gas, and that the wind height on the galaxy 
plane may still be related to the elapsed time, the result 
is that the metallicity of the X-ray emitting gas is tending 
towards lower abundances as one moves from the external 
to the inner regions. 

Two possible objections to this picture may be the con- 
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tradictory behaviour of S, and the fact that fragmenting may 
be not effective in increasing the clouds surfaces. E.g., if a 
spherical cloud fragments in N smaller spherical clouds and 
the density does not change, than the total surface S scales 
with S ocN 1/3 . 



6.3 Energetics of Charge-Exchange emission 

From a Suzaku observation, iTsuru et alj (|2007t ) calculated 
that the CE process can significantly contribute to the OVII 
A21-22 A lines observed from the region of extended emis- 
sion known as 'the Cap'. This emission is considered a part 
of the M82 outflow, and lies ~ 11' (11.6 kpc) north of the 
M82 nucleus. The best-fitting spectral models of the Cap 
give a density for the hot gas of ~ 5 • 10~ 3 cm -3 . Under the 
(extreme) hypothesis that all O ions suffer the CE process 
and emit O K lines, and assuming a density of the order 
10~ 3 cm" 3 for the neutral material, Tsuru et al. estimated 
an upper limit to the contribution from CE of 5 • 10~ 6 pho- 
tons s" 1 cm -2 , to be compared with an observed flux of 
6 • 1CT 6 photons s" 1 cm . To make a similar estimate for 
the M82 central region, one has to rescale for the volume 
and densities. The volume and the hot gas density are taken 
from Table El (1.4 kpc 3 and 32 • 10~ 3 cm" 3 ). The mass and 
density of the molecular material c ontained in the same vo l- 
ume are 3.7- 10 8 M© and 11 cm" 3 | Young fe Scovillell 19841 ). 
The rescaled upper limit on the expected flux would thus 
be 4.0 • 10~ 2 photons s" 1 cm" 2 , to be compared with the 
observed flux for the OVII triplet which is 2.4- 10~ 4 photons 
s" 1 cm" 2 . Thus the flux of the OVII triplet can be fully jus- 
tified if the CE process involves < 1% of the OVIII ions. Of 
course, this is still an upper limit, as the flux contributed to 
the triplet by thermal emission is highly uncertain. 



6.4 Reliability of Oxygen abundances. 

Discrepancies between measurements in stars 
and gas 

Taken at its face value, the X-ray derived O abundance is 
found to be smaller than the ones derived for Fe, Mg and 
other heavier a-elements. This effect is true across all the 
outflow, although it seems more severe in the central re- 
gions: using EPIC measurements, the O/Mg abundance ra- 
tio is 0.7 ± 0.3 in the external outflow, and 0.2 ±0.1 in the 
inner part. On the other hand, WW95 calculates a O/Mg 
abundance ratio of ~ 10-20 for all stars with masses between 
11 and 40 Mq. Similarly, the observed O/Si ranges from 0.2 
and 1.0 while WW9 5 has 2-10. This was already noticed by 
lUmeda et all (120021). who found in A SCA data 0/Mg~ 0.2 
and 0/Si~ 0.15. lUmeda et all (2002) had some success in 
explaining these abundance ratios by invoking hypernova 
explosions, i.e. energ etic (E > 10 52 erg) core-collapse su- 
pernovae with massive (M > 25Mq) metal-poor [Z = 0) 
progenitors. 

However, the X-ray derived oxygen absolute abundance 
is smaller than the IR-derived one. This looks like a more 
severe problem, because the IR-derived abundances should 
trace the composition of the gas before any enhancement by 
the last starburst occurs. 

A significant contribution from CE emission has been 
found in the O VII lines (Sect. l4~2)l . thus challenging also 



the O VIII lines because CE is expected to be f ound in 
lines from both ionisation states (|Diuric et al.ll2005T ). While 
this might question the reliability of O abundance measure- 
ments from X-ray observations, the contributions from CE 
and thermal emission should add together and result in a 
stronger line, and hence in a larger measured abundance. 
Thus it seems that the problem of O under-anbundance is 
onl y made worse. 

iMarcolini et~aT1 (|2005T ) found, in simulation of 
wind/cloud interaction, that the conductive interface 
between wind and clouds has a lower temperature than 
the wind bow shock. Thus they suggested that, since 
the emissivity of different ions is strongly temperature 
dependent, the X-ray line emission from low mass ions (O 
VII, O VIII) should be more influenced by cloud material, 
while the contribution from wind material should be larger 
for higher ionisation species (Mg and heavier). On top of 
this, one could also consider the CE as a cooling mechanism 
for O ions, because in reactions between O and an element 
X like 

O q+ + X -> O c,_1)+ * + X + (5) 

the O ion looses energy by falling in a lower ionisation state, 
and may eventually end out of the 'temperature window' 
available for emission in the 0.5-2.0 keV band. So the bulk 
of the O mass could probably be traced by measuring O in 
the UV and far UV bands. 



7 CONCLUSIONS 

We have reported on a deep (~ 100 ks) X-ray observation of 
the starburst galaxy M82 conduced with the XMM-Newton 
satellite. Spectra of both the galaxy centre and the outflow 
have been presented. 

By observing the M82 nuclear regions with EPIC, at 
least three spectral components were found: a power-law 
making most of the flux at energies larger than a few keV, a 
collisionally ionised plasma accounting for most of the flux 
at lower energies, and three narrow lines whose origin may 
be attributed to charge exchange emission. 

The power-law emission may be ascribed to the 
many point sources present in the M82 central region, 
among whom M82-X1 is by far the most luminous. Its 
data from this observa tion have been reported elsewhere 
jMucciarelli et al.ll2006l ). 

The thermal emission has been found to have a double- 
peaked differential emission measure, the first peak being at 
~ 0.5 keV and the second one at ~ 7 keV. A similar tem- 
perature structure, albeit with a smalle r temperature ratio , 
is predicted by the model proposed bv lMasai et al. I (|2002T ). 
According to these authors, if a plasma is subject to particle 
acceleration co-spatial with the thermal emission, then the 
energy spectrum of the electrons is not fully thermal but 
includes a non-thermal part, simulating a hot component in 
the X-ray photon spectrum. Indeed, the Fe line at ~ 6.7 keV 
is consistent with a plasma with temperature kT ~ 0.5 keV 
but not with a temperature of ~ 7 keV. 

The double-peaked structure of the DEM was confirmed 
in a series of spectra of the outflow taken with EPIC at in- 
creasing distance from the galactic plane. The plasma tem- 
peratures do not change significantly, the lower energy peak 
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moving at lower energies further from the galactic plane, 
while the upper energy peak stays constant across the cen- 
tre and all the outflow. 

Three lines were found in the EPIC spectra which may 
be attributed to charge exchange reactions involving neutral 
Mg and Si present in dust grains and the O VIII n = 3 — > 1 
transition. The equivalent width of these lines depend on 
the camera (MOS or pn) and on the spectral model used for 
the plasma (MEKAL or APEC), ranging from ~ 4 to ~ 40 
eV. The RGS data confirm the CE lines from Mg and Si, 
while the CE line from O is not confirmed. RGS data also 
hint for a CE origin of the O VII A21. 6-22.1 triplet (EPIC 
cannot place constraints on the line intensity ratios of this 
triplet because of insufficient spectral resolution). 

Both the O VIII line at 19 A and the O VII triplet 
at 21-22 A consist of two components: the main one, at a 
redshift similar to that of the other elements' lines, and a 
secondary one with lower intensity which exhibits a blueshift 
(z ~ —0.012 corresponding to Av ~ 3800 km s" 1 ). More- 
over, the lines making the dominant (redshifted) component 
O VII triplet have line intensity ratios consistent with CE 
emission but not with a thermal plasma. 

The chemical abundances of O (taken at its face value 
and ignoring any possible contribution from CE), Ne, Mg, 
Si and Fe as derived with EPIC increase when moving out 
of the galaxy centre and towards the external regions of 
the outflow, while the behaviour of S is unclear. On the 
other hand, since the RGS is a slitless spectrometer with a 
large field of view, the spectral properties derived with this 
instrument should be regarded as field-of-view-averaged. In- 
deed, the RGS-derived abundances are consistent with being 
a spatially-weighted average of the EPIC ones. Abundance 
ratios, such as Mg/Fe or Si/Fe also show spatial struct ure. 

Confirming our previous study l)Origlia et al.l[20o3 ). the 
X-ray derived O abundance in the M82 centre is lower than 
the IR-derived one, at variance with the expectation that 
the X-ray emitting gas would be enriched with respect to 
the stars. Reminding that the X-ray emission from different 
elements preferably occurs at different temperatures, and 
that the main site for X-ray emission is the i nterface between 
the h ot gas and the colder ambient material ( Marcolini ct al. 
120051 ). we suggest that a significant fraction of the O VIII 
and O VII ions already have cooled off through interaction 
with cold material or the CE mechanism. Observations in 
the far ultraviolet regime could help clarify this issue. 
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